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Abstract. Estimating black hole masses of blazars is still a big chal- 
lenge. Because of the contamination of jets, using the previously sug- 
gested size - continuum luminosity relation can overestimate the broad 
line region (BLR) size and black hole mass for radio-loud AGNs, in- 
cluding blazars. We propose a new relation between the BLR size 
and Hp emission line luminosity and present evidences for using it 
to get more accurate black hole masses of radio-loud AGNs. For ex- 
tremely radio-loud AGNs such as blazars with weak/absent emission 
lines, we suggest to use the fundamental plane relation of their ellipti- 
cal host galaxies to estimate the central velocity dispersions and black 
hole masses, if their velocity dispersions are not known but the host 
galaxies can be mapped. The black hole masses of some well-known 
blazars, such as OJ 287, AO 0235+164 and 3C 66B, are obtained using 
these two methods and the M - a relation. The implications of their 
black hole masses on other related studies are also discussed. 
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1. Introduction 

In the last two decades, dynamical measurements have clearly revealed that su- 
permassive black holes exist in the center of nearby galaxies. However, dynamical 
methods can not be applied to most of AGNs because they are too bright. Currently 
the most reliable method for AGN black hole mass estimation is the reverberation 
mapping. Using this technique, the BLR size can be measured using the time lag 
between the variations of continuum and emission line fluxes. The black hole mass 
can be derived from the BLR size and the characteristic velocity using the virial 
law. So far, reverberation mapping studies have yielded black hole masses of about 
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40 Seyfert 1 galaxies and nearby quasars ( Kaspi et al. 2000, 2005; Peterson 
et al. 2004). Using the observed data of these reverberation mapping AGNs, an 
empirical relation between the BLR size (R) and continuum luminosity at 5 100 A 
(L ° ) has been derived by Kaspi et al. (2000, 2005), which has been frequently 
adopted to estimate the BLR size and the black hole masses for large samples of 
AGNs, including radio-loud quasars. However, the optical continuum luminosity 
of radio-loud AGNs may not be a good indicator of ionizing luminosity. Powerful 
jets of blazar-type AGNs may significantly contribute to the optical continuum lu- 
minosity. Therefore, using the R — ^ 5100 ^ relation may significantly overestimate 
the actual BLR size and the black hole mass of radio-loud AGNs. In addition, an- 
other tight correlation between black hole mass and bulge velocity dispersion (a) 
has been found for nearby galaxies (Tremaine et al. 2002) and for a few Seyfert 
galaxies as well (Ferrarese et al. 2001). Such a relation suggests a possibility to 
estimate the black hole masses of AGNs using the measured bulge velocity dis- 
persions. Especially for BL Lacertae objects, the reverberation mapping technique 
cannot be applied because they show no or weak emission lines in their optical 
spectra. Using the Mbh-& relation may be the only way to derive their black hole 
masses, though measuring a is possible only for nearby sources. We have to seek 
for other methods for most blazars because the velocity dispersions of their host 
galaxies are too difficult to be measured. 

In this paper we report our progress on estimating the black hole masses of 
radio-loud AGNs using a new BLR size - Hp emission line luminosity relation 
and the fundamental plane relation of the elliptical host galaxies. The results on 
some well-known blazars are also presented. 



2. The BLR size - Hg luminosity relation 

Using the available data of BLR sizes and IL3 fluxes for 34 AGNs in the reverbera- 
tion mapping studies (Kaspi et al. 2000), we investigated the relation between the 
BLR size and the Hg emission line luminosity (Wu et al. 2004). An empirical 
relation between the BLR size and Hg luminosity was derived as: 

Log R (It - days) = (1.381 ±0.080) + (0.684±0.106)Log (L H/3 /10 42 ergss" 1 ). 

(1) 

The Spearman's rank correlation coefficient of this relation is 0.91. In Fig. 1 we 
show the dependence of the BLR size on Lh p and ^ 5100 ^- Obviously these two 
relations are similar, which means that the R — Lh^ relation can be an alternative of 
the R — L 5 Qo ° relation in estimating the BLR size for radio-quiet AGNs. We ap- 
plied both me°it — Lh» and R — L ° relations to estimate the black hole masses 
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of 87 radio-loud quasars and compare them in Fig. 2. Evidently the masses ob- 
tained with the R — Lh^ relation are systematically lower that those obtained with 
the R — L i ° relation for some extremely radio-loud quasars. The difference be- 
tween two hole mass estimates is smaller when the radio-loudness is small but 
becomes larger as the radio-loudness increases. For some individual quasars with 
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higher radio-loudness, the black hole mass estimated with the R — L ° relation 

5100A 

can be 3~10 times larger than that estimated with the R — Lu fi relation. Recently 
Kong et al. (2006) also extended such a study to the broad UV emission lines Mgll 
and CIV, and obtained the BLR size - UV emission line luminosity relations. Our 
results demonstrated that using the BLR size - emission line luminosity relations 
can avoid the overestimations of the black hole masses for blazar-like radio-loud 
AGNs. 
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Figure 1. The R — Lh p relation and the R — L ^ relation. The open and filled symbols 
denote Seyferts and quasars respectively. The figure is taken from Wu et al. (2004). 



3. Black hole masses estimated from the fundamental plane relation 
of AGN elliptical host galaxies 

For radio-loud AGNs such as BL Lacs with weak/absent emission lines, the emis- 
sion line based methods for black hole mass estimations can not apply. Directly 
measuring their stellar velocity dispersion is also difficult. However, the host 
galaxies of BL Lacs are virtually ellipticals. It is well known for ellipticals that 
three observables, namely the effective radius (R e ), the average surface brightness 
(< fi e >r in R-band) and the central velocity dispersion (a), follow a tight funda- 
mental plane relation. For about 300 normal ellipticals and radio galaxies, Bettoni 
et al. (2001) found that the fundamental plane can be robustly described as 

log R e = (1.27 ± 0.04) log a + (0.326 ± 0.007) < fi e > R -8.56 ± 0.06, (2) 

This relation provides us another way to estimate the central velocity dispersions 
and then the black hole masses of AGNs (Wu, Liu & Zhang 2002). 
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Figure 2. Comparison of the black hole masses of radio-loud quasars estimated with two 
R-L relations, and the dependence of the black hole mass difference on radio loudness. 
The figure is taken from Wu et al. (2004). 
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Figure 3. Histograms of the derived black hole mass distribution of HBLs and LBLs. The 
figure is taken from Wu et al. (2002). 

Using the imaging data of BL Lacs obtained from the HST snapshot survey 
(Urry et al. 2000), we adopted the fundamental plane relation to estimate the 
central velocity dispersions and black hole masses of 51 high-frequency peaked 
BL Lacs (HBLs) and 12 low-frequency peaked BL Lacs (LBLs). Our results show 
no significant difference in the black hole masses between HBLs and LBLs (see 
Fig. 3). We also applied the same method to 10 radio galaxies (RGs), 10 radio- 




Figure 4. Histograms of the derived black hole mass distribution of radio galaxies, radio- 
loud and radio-quiet quasars. The figure is taken from Wu et al. (2002). 



loud quasars (RLQs) and 13 radio-quiet quasars (RQQs) which have been imaged 
by HST (McLure et al. 1999), we found that there are no significant differences 
in the black hole masses among these different types of AGNs with elliptical host 
galaxies (see Fig. 4). 



4. Black hole masses of three well-known blazars 

4.1 AO 0235+164 

Liu, Zhao, & Wu (2006) applied the R-L Hfj relation suggested by Wu et al. (2004) 
to a well studied BL Lac object AO 0235+164 and estimated its black hole mass 
as 5.8 x 10 8 M Q , which is consistent with the mass (3.6 x 10 8 M©) estimated from 
the Mbh — & relation by taking the narrow emission line width as a surrogate of 
a, but is much smaller than the mass (1.5 x 10 9 M Q ) obtained from the R — L ^ 
relation. This again demonstrates that using the previous R — ^ relation can 
overestimate the black hole masses of blazar-like AGNs. Our stud°y indicates that 
the black hole mass of AO 0235+164 is most likely around 5 x 1O 8 M . 



4.2 OJ287 

We applied the fundamental plane relation of ellipticals to a well-known BL Lac 
object OJ 287, which is a possible supermassive black hole binary system. Liu & 
Wu (2002) estimated its primary black hole mass to be about 4 x 10 8 M©, which is 
consistent with the upper limit (10 9 M©) obtained by Valtaoja et al. (2000) based 
on a binary black hole model for OJ 287. This value is also within the range of 
typical black hole masses for BL lac objects (Wu, Liu & Zhang 2002). We noticed 
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that our derived black hole mass is much smaller than 1.8 x 10 M , which is 
required in a new precessing binary black hole model of OJ 287 (Valtonen 2007). 

4.3 3C 66B 

3C 66B is a well-known nearby FR I radio galaxy. The central velocity dispersion 
of 3C 66B has been obtained by Balcells et al. (1995) as 348±29 km/s based on 
their spectroscopic observations. Using the M - a relation suggested by Tremaine 
et al. (2002), we estimate that log(M/M ) = 9.097 ± 0.175, namely the mass 
of the central black hole of 3C 66B is about 1.25 x 10 9 M o . This is within the 
range of typical black hole masses for radio galaxies (Wu, Liu & Zhang 2002), 
but one-order of magnitude smaller than the maximum mass estimated by Sudou 
et al. (2003), 5 x 10 10 M Q , which was obtained by modeling the observed orbital 
motion with a binary black hole model. A very recent study by Iguchi et al. (2010) 
supported that the larger black hole mass is (1.21q 2) x 10 9 ^/©> which is consistent 
with the value obtained by the M - a relation. Janet et al. (2004) adopted the black 
hole mass of 3C 66B obtained by Sudou et al. (2003) to calculate the gravitational 
wave emission from 3C 66B and argued that the binary black hole model can be 
excluded because observations from the Pulsar Timing Array (PTA) did not detect 
the gravitational wave from 3C 66B. However, if the smaller black hole mass is 
adpted, the gravitational wave emitted from 3C 66B will be weaker than the current 
PTA detection limit and the binary black hole model can not be excluded. 

5. Summary 

We proposed to use the BLR size - emission line luminosity relation and the fun- 
damental plane relation of the elliptical host galaxies to estimate the black hole 
masses of radio-loud AGNs. We demonstrated that with the first relation we can 
get more accurate black hole mass estimates for radio-loud AGNs than using the 

usual R — L ° relation, and for some radio-loud AGNs such as BL Lacs the sec- 

5100A 

ond method is probably the only available one for their black hole mass estimations 
in the case when directly measuring the stellar velocity dispersions is difficult. We 
have adopted these methods, as well as the M - a relation, to estimate the black 
hole masses of three well-known blazars, AO 0235+164, OJ 297 and 3C 66B, and 
the results are helpful to other related study about their physics nature. Finally, 
we would like to mention that these two methods can be also applied to estimate 
the black hole masses of high redshift AGNs, including blazars, with high quality 
spectroscopic and imaging observations. 

This work is supported by the NSFC Grants (No. 10473001, 10573001, 10525313 
and 11033001) and the 973 program (No. 2007CB8 15405). 
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